Trembling motion of relativistic electrons in a magnetic field 
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Zitterbewegung (ZB, the trembling motion) of free relativistic electrons in a vacuum in the pres- 
ence of an external magnetic field is calculated. It is shown that the motion of an electron wave packet 
has intraband frequency components, corresponding to the classical cyclotron motion, and several 
interband frequency components corresponding to the Zitterbewegung. For a two-dimensional situ- 
ation, the presence of a magnetic field makes the ZB motion stationary, i.e. not decaying in time. 
We show how to simulate the ZB in a magnetic field using trapped ions and laser excitations in 
the spirit of recently observed proof-of-principle ZB simulation by Gerritsma et al. Nature 463, 68 
(2010). It is demonstrated that, for the parameters of the Dirac equation simulated by the above 
experiment, the effect of a magnetic field on the Zitterbewegung is considerable. 

PACS numbers: 31.30.J-, 03.65.Pm, 41.20.-q 



The phenomenon of Zitterbewegung was theoretically 
devised by Schrodinger [l| for free relativistic electrons in 
a vacuum. Schrodinger showed that, due to a noncom- 
mutativity of the quantum velocity operators with the 
Dirac Hamiltonian, relativistic electrons experience the 
trembling motion even in absence of external potentials. 
Here one deals with a purely quantum effect since it goes 
beyond Newton's first law of classical motion. It was 
later recognized that the phenomenon of ZB is due to an 
interference of electron states with positive and negative 
electron energies 0, Q . The predicted frequency of ZB 
oscillations is very high, corresponding to hujz — 2mc^, 
and its amplitude is very small being around the Comp- 
ton wavelength h/mc = 3.86x 10"'^ A. Thus, it is im- 
possible to observe this effect using present experimental 
arrangements. In fact, even the principal observability of 
ZB was often questioned in the literature, see e.g. 0, [1]. 
However, in a very recent paper Gerritsma et al. Q sim- 
ulated the 1-f 1 Dirac equation (DE) and the resulting 
Zitterbewegung with the use of trapped Ca ions excited 
by appropriate laser beams. The remarkable advantage 
of this method is that one can simulate the basic param- 
eters of DE, i.e. mc^ and c, and give them desired values. 
This results in a much lower ZB frequency and a much 
larger ZB amplitude. These were in fact observed. 

The purpose of our work is twofold. First, we calculate 
the Zitterbewegung of relativistic electrons in a vacuum 
in the presence of an external magnetic field. The pres- 
ence of a constant magnetic field does not cause electron 
transitions between negative and positive electron ener- 
gies. On the other hand, it quantizes the energy spec- 
trum into Landau levels which brings qualitatively new 
features into the ZB. The problem of ZB in a magnetic 
field was treated before ^ , but the results were limited to 
the operator level and suffered from various deficiencies 
which we briefly mention below. A similar problem was 
treated in Ref. [8| on the operator level in weak magnetic 
field limit. Our second purpose is to show how to simu- 
late the Dirac equation for electrons in a magnetic field 
with the use of trapped ions. We demonstrate that, using 



the parameters of DE simulated in Ref. 6] , the effects of 
the simulated magnetic field on ZB are sufficiently strong 
to be observable. 

The Hamiltonian for a relativistic electron in a mag- 
netic field is 



(1) 



where -k = p — qA is the generalized momentum, q is 
the electron charge, ai and f3 are the Dirac matrices in 
the standard notation. Taking a magnetic field B\\z we 
choose the vector potential A = {—By, 0, 0) and look for 
eigenfunctions in the form 



ijf(r) = e*'=-^+*'^-^$(y). 



(2) 



Introducing the magnetic radius L = y^h/eB and ^ = 
y/L — kxL, we have y = S^L -\- k^L^, eBy/h = y/L"^, and 
d/dy = {l/L)d/d£,. Defining the standard raising and 
lowering operators for the harmonic oscillator = (^ -I- 
d/dO/V2 and a+ = (^ - d/d^)/V2 one has [ay, a+] = 1 
and ^ = {ay + a+)/\/2. Now the Hamiltonian reads 



n = 



where 1 is the 2x2 identity matrix, and 
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with uj = ypicjL. An eigenstate |n) of H is character- 
ized by five quantum numbers: n, kx,kz, e, s, where n is 
the Landau level number, k^ and kz are the wave vec- 
tor components, e = ±1 labels the positive and negative 
energy branches, and s = ±1 is the spin index. In the 
representation of Johnson and Lippmann Q the state |n) 
is 



Su{£Enk, + mc^) \n - 1) 

siieEnk, + mc^) \n) 

[suhkzC- sihhJn) |n - 1) 

-{SuhUn + sihkzc) In) 



, (5) 
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where Su = {s + l)/2 and s; = (s — l)/2 are the projec- 
tion operators on the states s = ±1, respectively. The 
frequency is a;„ = LUy/n, the energy is 



nkz 



(6) 



and the norm is N^ek^ — C^E^k + 2emc^i?„fc^ )^/^. In 
this representation the energy Enk^ does not depend ex- 
phcitly on s. The harmonic oscillator states {r\n) have 
the standard form. 

In order to calculate the average electron posi- 
tion we introduce four-component operators A ~ 



diag(ay, 



ay,ay) and = diag(a+, a+, a+ 



a+). Now 



we define y = L{A + A+)/V2 and X = L{A - A+)/iV2 
in analogy to the position operators y and x. We use 
averaging of A and A'^ operators in the Heisenberg pic- 
ture A{t) = e'"^/^A{Q)e-^"^/^ over a wave packet fir). 
Inserting two unity operators 1 — |ii)(n| one obtains 



{A{t)) =y{,f\n)e 



ieEnt/h. 



nn' 



n|i|n') 



-ie'E. 



"'*/^n'|/), (7) 



and similarly for {A'^ (t)). The wave packet is assumed to 
be in the form /(r)* — fl{z)f^y{x,y), where i = 1,2,3,4 
labels the packet components. Below we limit our cal- 
culations to a packet with the second nonzero compo- 
nent. The summation over n, n' denotes summations over 
n,n',e,e',s,s' and integrations over ^ dk^ ^ dk^ , dk ^ . 
We obtain explicitly 

/oo 
dk^dk'MkMik'.) X 

Ei{tE„k^-t' E^,^, )t/h 
e " ^ Xnek^Xn' e k'^ ^ 



oo 



E 



dk,dk,sis[Fl*{k,)Fl,{kMA\r^), (8) 



where Xnek^ = [e-E^k, + ■m<?)l'Nnek,, and 
in which 
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5L(fcx,y)e-^« H„(Ody, (9) 
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fUx,y)e^'-^dx, (10) 



and 



n{z)e'''''dz. 



(11) 



The selection rules for (n|^|n') are n' — n + 1, k'^ ^ k^, 
k'^ = kz, while for (n|^+|n') they are n' — n~\, k'^ — k^, 
k'^ — kz- After some manipulation we finally obtain 

(Mt)) - ^ E {It + ic + lit + lis) , 

n 

{A+{t)) =\Y1 ^f^^^ul,,n {It + Ic lit ll7) , 
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FIG. 1: Calculated motion of the electron wave packet in vac- 
uum in the presence of a gigantic magnetic field B = 2x lO" T 
according to 3-f 1 Dirac equation for two packet widths. Both 
interband and intraband frequencies are present in the spec- 
trum. Time is given in tc = h/mc^, positions are given in 
Compton wavelengths Ac = h/mc. 
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and 



^ m.n 



F::{k,)Fi{k^)dk^, 



(16) 



for i,j — 1,2,3,4. We perform specific calculations tak- 
ing f{r) in the form of an ellipsoidal Gaussian packet 
(0, /(r), 0, 0) characterized by three widths dx, dy, dz and 
having a non-zero momentum hko — ?i(fcox, 0, 0). The 
quantities gxy{kx,y), Fl{kx), and U^^^ can be obtained 
analytically, see Ref. [l(]| . 

In the non-relativistic limit: hiu^cpz ^ mc^, equa- 
tions ([T^ — reduce to the cyclotron motion with 
the frequency Wc = {En+i,k, — En,k,)/fi — eB/m and the 
radius r = ko^L^ = mv/eB. In this limit the ZB part 
of the motion has the amplitude several orders of mag- 
nitude smaller than Ac. In the opposite relativistic limit, 
in which ftw, cpz — rac^ , both the cyclotron and ZB com- 
ponents have comparable amplitudes. In Fig. [1] we plot 
the average position in x and y directions for an electron 
wave packet in a gigantic magnetic field B = 2x10^ T 
for two packet parameters. It is seen that the ZB oscil- 
lations consist of several frequencies. This is the main 
(12)efrect of the magnetic field, which quantizes both posi- 
tive and negative electron energies into Landau levels. In 
^3^1arger time scale the oscillations in 3-f 1 space go through 
decays and revivals, but finally disappear. 
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The involved DE must have at least 2+1 dimensions 
since a magnetic field B \\ z couples the electron mo- 
tion in X and y directions. The above equations, derived 
for the 3+1 DE, can be reduced to the 2+1 case by set- 
ting \gj{kz)\'^ — S{kz). The main difference between 3+1 
and 2+1 cases is that for the 3+1 case the ZB of the 
wave packet has a transient character (it decays in time) , 
whereas in the 2+1 case it has a permanent character go- 
ing through decays and revivals. This difference, due to 
the integration over in the 3+1 case, is in agreement 
with general predictions of Lock [llj. However, a slow 
decay of ZB in time will occur also in the 2+1 case since 
the trembling electron will emit radiation. This is possi- 
ble because a Gaussian wave packet is not an eigenstate 
of the Dirac Hamiltonian. Also a broadening of Landau 
levels due to external perturbations results in a transient 
character of ZB, c.f. |12| . 

The main experimental problem in investigating the 
above ZB phenomenon in an external magnetic field is 
the fact that for free relativistic electrons the basic ZB 
(interband) frequency corresponds to the energy hujz — 
1 MeV, whereas the cyclotron frequency for a magnetic 
field of 100 T is hojc — 0.1 eV, so that the magnetic effects 
in ZB are very small. However, as mentioned above, it is 
possible by now to simulate the Dirac equation changing 
thereby its basic parameters. This gives a possibility to 
modify drastically the ratio huJc/^mc^ making its value 
much more advantageous. To simulate the DE for an 
electron in a magnetic field, we transform it first to the 
form Ti, = c otiPi + 5mc^ , using the unitary operator 
P = (5((5 + (3)/\/2, where 6 = axUyazP jl3]. After the 

H 
^ ij't 

where 



transformation the Hamiltonian is "H = 



H = 



CPz 
CPx 



CPx 
-CPz 



hiua,, 



(17) 



We then follow the approach of Lamata et al. [ij] and 
consider a four-level system of Ca or Mg ions. Simu- 
lations of cpx and cpz terms in the above Hamiltonian 
are realized the same way as for a free Dirac parti- 
cle. Thus one uses pairs of the Jaynes-Cumminngs (JC) 
interaction: -ffj^ = hiriVt{a^ ae^'^'" + o'^a+e^*'^'^), and 
the anti- Jaynes-Cumminngs (AJC) interaction: H'j^j^j — 
hr]fl{a~^ a'^ e'^'^'' + a~ ae~'^'^''), while a simulation of mc^ 
is done by the carrier interaction H'^ — hVl{a^e'^'^'' + 
(T~e~"^'=). Here Q. and are the coupling strengths and 
is the Lamb-Dicke parameter. On the other hand, a sim- 
ulation of CLy and (which include the magnetic field) 
can be done by single JC or AJC interactions. Within 
the notation of Refs. [T3 - [T6t one may simulate H by the 
following set of excitations 



a^{bc) 



H 



JC{ad) 
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^ ^AJC(bc) 



where H^] = 2T]qQajAqPq, pq = ih{a+ - aq)/ Ag, j,q = 
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FIG. 2: Calculated motion of the electron wave packet as 
described by 2+1 DE, simulated by trapped ^"Ca^ ions for 
three values of the effective rest energy Q. Trap parameters: 
rj = 0.06, = 27r X 68 kHz, A ~ 96A, fcox = A"\ dy = A^/2, 
dx = 0.9dy. Simulations correspond to ratios k — hujc/2mc^ = 
16.65(a), 1.05 (b), 0.116 (c), respectively. Positions are given 
in L = V2A. 



X, z, and the Pauli matrices aj are obtained from com- 
binations of Hj'fj and H'^j^j with appropriate phases 
<j)r and (/){,. Here the spread of the ground ion wave 
function is = ^hjlMvq and the Lamb-Dicke pa- 
rameter is r\q — k^hjlMvq^ where M is ion's mass, 
Vq is trap's frequency in the q direction and fc is the 
wave vector of the driving field in a trap. The sub- 
scripts in parenthesis of Eq. (IT5|) denote states involved 
in a given transition. The JC interaction gives hy in 

11^2 and fly in li2\ elements of the Hamiltonian H in 

Eq. (fTT]) . while AJC gives hy in a+ in H^^ el- 

ements. The simulation of 3+1 DE by Eq. (|18l) can be 
realized with 12 pairs of laser excitations. If one omits 
the Pz interaction, which corresponds to the 2+1 DE, one 
needs 8 pairs of laser excitations. The simulated mag- 
netic field can be found from the correspondence (see 
Ref. [l|): fly - a+ = y/2L{d/dy) = 2A{d/dy), which 
gives A <!=> L/^/2. Since c <^ 2r/Afi and mc^ hfl, we 
have 



heB _ -qn 
m{2mc^) \ Vl 



(19) 



Therefore, by adjusting frequencies VL and £7 one can 
simulate different regimes of the ratio k — hui^j^mc? . 
In Fig. [5] we show the calculated ZB for three values of 
k: 16.65, 1.05, 0.116. It is seen that, as k gets larger 
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(i.e. the field intensity increases or the effective gap de- 
creases), the frequency spectrum of ZB becomes richer. 
This means that more interband and intraband frequency 
components contribute to the spectrum. Both intraband 
and interband frequencies correspond to the selection 
rules n' = rt ± 1 so that, for example, one deals with 
ZB (interband) energies between Landau levels n = 
to n' = 1, and n = 1 to n' = 0, as the strongest con- 
tributions. For high magnetic fields the interband and 
intraband components are comparable, so one can legiti- 
mately talk about the ZB. We believe that ZB oscillations 
of the type shown in Fig. [2^, based on the 2-1-1 DE for 
K — fujjc/2mc? ^ 1, are the best candidate for an ob- 
servation of the simulated trembling motion in the pres- 
ence of a magnetic field. The calculated spectra shown in 
Fig. [5] use the simulated parameters already realized ex- 
perimentally, see @. Notice the tremendous differences 
of the time and position scales between the results for free 
electrons in a vacuum shown in Fig. [T] and the simulated 
ones in Fig. [5] 

The anisotropy of ZB with respect to the {x{t)) and 
{y{t)) components, seen in Figs. [Hand [21 is due to the 
initial conditions, namely k^x 7^ and k^y = 0. A simi- 
lar anisot ropy was predicted in the zero-gap situation in 
graphene (lOj . 

We emphasize that the Dirac equation llj and our re- 
sulting calculation, as well as Eq. ([TT]) and its simula- 



tion in Eq. ([18]), represent the 'empty' Dirac Hamilto- 
nian which does not take into account the 'Fermi sea' 
of electrons with negative energies in a vacuum. This 
one-electron model corresponds to the original consider- 
ations of Schrodinger's. On the other hand, the filled 
states of electrons with negative energies may affect the 
phenomenon of ZB, see (sl. Il7l|. 

Finally, one should recognize that the experiment of 
Gerritsma et al. @ simulates not only the 1-t-l Dirac 
equation for free relativistic electrons in a vacuum but 
also the two-band k.p model for electrons in narrow-gap 
semiconductors and the ZB resulting from this descrip- 
tion [3, In fact, the results of Ref. [1] look remark- 
ably similar to our predictions [l^. 

In summary, we calculated the trembling motion of rel- 
ativistic electrons in a vacuum in the presence of a mag- 
netic field for 3+1 and 2+1 spaces. In contrast to the 
no-field case, the presence of a magnetic field results in 
many interband frequencies contributing to the trembling 
motion. In the 2+1 case the ZB oscillations of the elec- 
tron wave packet are stationary, i.e. they do no decay in 
time. We indicate how to simulate the Dirac electron in 
a magnetic field and the resulting ZB using trapped ions 
and laser excitations. We show that, for the parameters 
of DE simulated very recently by Gerritsma et al. Q , the 
effect of a magnetic field on the trembling motion should 
be clearly observable. 



[1] E. Schrodinger, Sitzungsber. Preuss. Akad. Wiss. Phys. [9 
Math. Kl. 24, 418 (1930). Schrodinger's derivation is re- 
produced in A. O. Barut and A. J. Bracken, Phys. Rev. [10 
D 23, 2454 (1981). 

[2] J. D. Bjorken and S. D. Drell, Relativistic Quantum Me- [11 
chamcs (McGraw-Hill, New York, 1964). [12 

[3] B. Thaller, The Dirac Equation (Springer- Verlag, Berlin, 

1992). [13 

[4] K. Huang, Am. J. Phys. 20, 479 (1952). 

[5] P. Krekora, Q. Su, and R. Grobe, Phys Rev. Lett. 93, [14 
043004 (2004). 

[6] R. Gerritsma, G. Kirchmair, F. Zahringer, E. Solano, R. [15 
Blatt and C. F. Roos, Nature 463, 68 (2010). 

[7] A. O. Barut and W. D. Thacker, Phys. Rev. D 31, 2076 [16 
(1985). This calculation of ZB in a magnetic field is lim- 
ited to operators, only four interband and intraband fre- [17 
quencies are predicted and the results diverge at mag- 
netic fields above the critical field heB/m > mc? , which [18 
corresponds to Bcr ^ 4: Ax 10^ T. [19 

[8] M. Villavicencio and J. A. E. Roa-Neri, Eur. J. Phys. 21, 
119 (2000). 



M. H. Johnson and B. A. Lippmann, Phys. Rev. 76, 828 
(1949). 

T. M. Rusin and W. Zawadzki, Phys. Rev. B 78, 125419 
(2008). 

J. A. Lock, Am. J. Phys. 47, 797 (1979). 

T. M. Rusin and W. Zawadzki, Phys. Rev. B 80, 045416, 

2009. 

R. E. Moss, and A. Okninski, Phys. Rev. D 14, 3358 
(1976). 

L. Lamata, J. Leon, T. Schatz, and E. Solano, Phys. Rev. 
Lett. 98, 253005 (2007). 

M. Johanning, A. F. Varron, and C. Wunderlich, J. Phys. 
B 42, 154009 (2009). 

D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. 
Mod. Phys. 75, 281 (2003). 

A. O. Barut and S. Mahn, Rev. Mod. Phys. 40, 632 
(1968). 

W. Zawadzki, Phys. Rev. B 72, 085217 (2005). 

W. Zawadzki and T. M. Rusin, arXiv:cond- 

mat/0909.0463 (2009). 



